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We study the interaction between two B — 1 states in the Chiral-Dilaton Model with scale 
invariance where baryons are described as non-topological sohtons arising from the interaction of 
chiral mesons and quarks. By using the hedgehog solution for the B = 1 states we construct, via 
a product ansatz, three possible B — 2 configurations to analyse the role of the relative orientation 
of the hedgehog quills in the dynamics. We investigate the behaviour of these solutions in the 
range of long and intermediate distances between the two solitons. Since the product ansatz breaks 
down as the two solitons get close, we explore the short range distances regime by building up a six 
quarks bag and by evaluating the interaction energy as a function of the inter-soliton separation. 
, We calculate the interaction energy as a function of the inter-soliton distance for the B = 2 system 

■ and we show that for small separations the six quarks bag, assuming a hedgehog structure, provides 

a stable bound state that at large separations connects with a special configuration coming from 
^ I the product ansatz. 
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I. INTRODUCTION 



The study of properties of nuclear matter under extreme conditions, at high temperature in relativistic heavy-ion 
physics and/or at high density as in compact stars, is one of the most challenging issues in hadronic physics. At 
the moment the phase diagram of hadronic matter is far away from being completely defined; in fact the regime of 
O I moderate temperatures and high densities, where the non-perturbative aspects of QCD dominate, seems to present 
^ ■ a rich set of scenarios and phases [ij. Because of the well-known "sign problem" in Lattice QCD one cannot 
" I , explore the whole phase space by the direct calculations from the fundamental theory with quark and gluon degrees 
of freedom. Thus, we must rely on effective field theories which are defined in terms of hadronic fields Q and that 
~^ ,, include the degrees of freedom and the symmetries relevant to the energy scale under study. In particular, chiral 
K*" ■ Lagrangians have been used in the past to describe the properties of the nucleons 0-01 and these models succeed 
in providing a reasonable description of the hadronic properties in the vacuum. On the other hand, the problem 
of studying nuclear matter at finite density by using these chiral Lagrangians does not lead to successful results. 
The failure of models based for instance on the linear cr-model is due to the constraints on the scalar field dynamics 
• imposed by the mexican hat potential [7l|. The solution to this problem, within models with a linear realization of 
^ ^ chiral symmetry, is to modify the meson potential, by incorporating another fundamental symmetry of QCD, namely 
, scale invariance [8l-[ll|. It has been already shown that with hadronic degrees of freedom, the Chiral-Dilaton Model 
' (CDM) provides a good description of nuclear ph ysic s at densities around po and it describes the gradual restoration 
• • ' of chiral symmetry at higher densities [T2j . In |l3l | the authors used this model interpreting the fermions not as 
nucleons, but as quarks and applying the Wigner-Seitz approach [l^ . showed that the new potential allows to reach 
densities higher than the ones obtained in the linear cr-model. 
• The Wigner-Seitz approximation [3], coming from solid state physics, has been used as a useful tool to mimic a 
dense system for both non-topological soliton models flsl-flTj, chiral soliton models (T7l - [20j and Skyrme models [2TI.[2^. 
Specifically the Wigner-Seitz method consists in replacing a cubic lattice by a spherical symmetric one where each 
soliton sits on a spherical cell of radius R with specific boundary conditions imposed on fields at the surface of the 
sphere. The configuration of the meson fields, centered at each lattice point, generates a periodic potential in which 
the fermions move. The Wigner-Seitz lattice represents a group of static field configurations with a precise symmetry 
connected to the boundary conditions imposed on the edge of each cell. Since the soliton solution contains also isospin 
degrees of freedom, it would be interesting to check if rotating the relative isospin orientation of the solitons we obtain 
a lower energy configuration, a study for which the Wigner-Seitz approximation is not well suited. 

In this work we present a first attempt to go beyond the Wigner-Seitz approximation and to study the behaviour 
of the baryon number B = 2 system using the CDM. We study the interaction energy of the system, defined as 
Vint — Eb=2 — '2.Eb=i, as a function of the distance d between the two solitons and discuss its behavior looking 
for the possibility of finding bound states. To develop this analysis we proceed by constructing a. B = 2 solution 
which describes well the long range behavior of the interaction. At the same time we build another solution which 
instead describes the intermediate/short range behavior of the interaction. By matching these two solutions at the 
intermediate range we construct the B = 2 interaction valid over the whole range. 

When the two solitons are very far apart we use the product ansatz and we study the interaction of the two solitons 
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by changing the relative orientation of the hedge hog qu ills. The product ansatz approach has been used in the past 
in the Skyrme model to describe the deuteron [23l - t26| . All these works start from the hedgehog solution in the 
5 = 1 sector and then the product ansatz is used to build different diskyrmions configurations where the isospin 
configuration plays a crucial role in providing the lowest energy state. The results given by the Skyrme model in this 
scheme are quite satisfactory and lead to a suitable description, concerning the quantum numbers, of the deuteron as 
a skyrmion. Nevertheless, it has been shown that using a more elaborate ansatz in which the two skyrmions merge 
in a toroidal configuration, the minimal energy state has axial symmetry [25l - l27| . More sophisticated approaches, 
including specific symmetries of the lattice, have been adopted to build skyrmionic matter and to describe on one side 
multi-baryon systems [1^ and on the other to study hot and dense baryonic matter [H, [s^l • 

The main difference of our approach with respect to the latter is the inclusion in the model of quark degrees 
of freedom which produce a repulsion at short distances associated to the Pauli principle. In order to overcome 
this problem here we build a six quarks bag where quarks are free to fill also the excited states. In soliton models 
including quark degrees of freedom, the B = 2 system has been studied starting from the self-consistent solution of 
field equations for six quarks, both assuming a hedgehog structure [3l|, Ull and a more sophisticated ansatz 33 1. In 
these works [H, |3^ the authors focused on a semi-bosonized version (similar to NJL models) of the chiral soliton 
model and they searched for topological solutions. The self-consistent solution of the field equations assuming a 
hedgehog configuration does not lead to a bound state in the B = 2 system and to avoid this problem the analysis is 
extended to an axially symmetric configuration of the meson fields leading to a stable system. 

In here instead we are looking for non-topological solutions of a CDM which explicitly includes quarks Hence, 
first of all we evaluate the energy for the CDM with six quarks at zero separation; the exact solution provides a stable 
solution with Vint ~ —8 MeV. Next we apply a variational calculation to the total energy of the bag in order to obtain 
an energy functional which varies with the distance d. This approach permits to join the results obtained in the long 
distances regime and the ones coming from this short distance calculation. 

The structure of the paper is the following. In Sec.|lT]we describe the model we are using, in Sec. Illll we present the 
calculations for the B = 2 system. In particular in Sec. IIII Al we show the long distances calculation method with the 
product ansatz, while in Sec. IIII 5] we present the results for the six quarks bag in the short range regime. Finally, in 
Sec. II VI we present our conclusions and future outlooks. 



II. THE MODEL 



The Lagrangian we are using, besides chiral symmetry, contains another fundamental symmetry in QCD, namely 
scale invariance. This invariance under dilatation is spontaneously broken in QCD because of the presence of the 
parameter Aqcd coming from the renormalization process. Formally, the non conservation of the dilatation current 
is directly connected to a not vanishing gluon condensate: 

{d,J^CD) = ^{F;Ax)F^^'^{x)). (1) 

In the papers of Schechter [33|, and Migdal and Shifman [slj they introduced a scalar field representing the gluon 
condensate, and whose dynamics is regulated by a potential chosen so that it reproduces (at mean-field level) the 
divergence of the scale current that in QCD is given by Eq. ([T]). The potential of the dilaton field is hence determined 
by the equation: 

^,^-4n^)-<^f =4e„..(i-)' (2) 

where the parameter e„ac represents the vacuum energy. The inclusion of massless quarks leads to a generalization, 
proposed in Ref. [s^, where also chiral fields contribute to the trace anomaly. In this way the single scalar field of 
Eq. ^ is replaced by a set of scalar andpseudoscalar fields {cr, tt, 0}. 

In our analysis we use, following Refs. [8l-ll2l|. a simplified chiral Lagrangian which contains only chiral fields coupled 
to quarks: 

C = ~ g^ia + in ■ rj^M + ^{d^crd'^a + d^n ■ d^n) - ^(00, a, tt). (3) 

The potential is given by: 



Via, tt) = Xlia^ + tt^) - ln{a^ + n^) - a^mla (4) 
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where: 



Xi-l^^^ = limi+mi) (5) 



Xl = \BS4 = ^{ml-ml) (6) 
ei=mlal (7) 

Here a is the scalar-isoscalar field, tt is the pseudoscalar-isotriplet meson field, (p is the dilaton field which, in the 
present calculation, is kept frozen at its vacuum value (/)o and ip describes the isodoublet quark fields. The decision to 
keep the dilaton field frozen is based on the results obtained in [l^] , where it has been shown that at low temper- 
atures the dilaton remains close to its vacuum value even at large densities. The vacuum state is chosen, as usually, 
at cr^ = and tt = 0. In the following calculations we use: — 5, mo- — 550 MeV, m^r = 139 MeV and /tt = 93 MeV. 

The self-consistent B = 1 solution for the Lagrangian in eq. ([3]) has been obtained by adopting the hedgehog ansatz, 
namely: 

crB=i(r) = ah{r) , 7rB=i(r) = h{r)f , 

^s=i(r) = ^f. A^i\ui)-\dt)). (8) 

The field equations for the Dirac components become: 

du 

— = g^hu+ (-€- gj,ah)v , (9) 

2 , , , 

— = — V - gr,hv + [e ~ g-^ahju , (10) 
ar r 

while for the chiral fields, the equations become: 

dr^ r dr 47r 9(T/i 

fh = -l^ + lh+^uv + — (12) 
dr-^ r dr 27r 9/i 

The solution has been obtained by imposing the following boundary conditions at r = 0: 

m'(0) = v{0) = 0, 

a^O) = h{0) = 0, (13) 

and at infinity: 

ah{oo) = /i(oo) = 0, 



where e is the quark eigenvalue. 

The energy of the soliton in the hedgehog configuration reads: 



where the quark-mesons interaction and the quark kinetic energies are: 



(14) 



£^5=1 = 47r / r^dr{E,nt{r) + Sfc™,Q(r) + E,{r) + E^{r) + Epot{r)) (15) 
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Eint{r) = — [g^f^hiu^ - v^) + 2g^huv] (16) 
^ , , 3 / dw d-u 2 \ 

Ekin,Q[r) = — [u- V— + -UV] (17) 

47r \ dr dr r J 
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and the energy density of the meson fields and of the potential read: 

2 



doh 
dr 

dr 



Epot{r) = V{(l)(i,(7h,h) 



(18) 

(19) 
(20) 



In Table H] we show all the contributions to the total energy of the soliton, coming from mesons and quarks, for the 
set of parameters: rric — 550 MeV, ra-^ — 139 MeV and g — 5. 



Table I. Contributions to tlie soliton total energy at mean-field level. All quantities are in MeV. 



Quantity 


Our Model 


Quark eigenvalue 


83.1 


Quark kinetic energy 


1138.0 


Ecr (mass+kin.) 


334.5 


Ett (mass+kin.) 


486.0 


Potential energy a — n 


88.0 


Eqa 


-101.4 


EqTT 


-787.0 


Total energy 


1157.8 



III. THE B = 2 SYSTEM 



To study the soliton-soliton interaction as a function of the distance d between solitons we proceed to construct 
a system with baryon number B ^ 2. At the moment, within the CDM, the finity density analysis of the chiral 
Lagrangian in eq.([3|) has been developed only using the Wigner-Seitz approximation in Ref. [l^, where the authors 
also included the vector mesons. The Wigner-Seitz approach, besides representing a first step in the study of dense 
matter, nevertheless implies the analysis of a lattice configuration of static solitons with a precise symmetry imposed 
by the boundary conditions on the surface of the spherical cell. 

To study the behaviour of the energy functional as a function of the separation distance between solitons, we proceed 
as follows. We construct a product ansatz solution valid for long distances; next we construct a. B = 2 hedgehog 
solution which holds at short distances. Finally we match the two solutions in the intermediate region by keeping 
that with minimal energy. 



A. Long and intermediate distances: the product ansatz approach 



The hedgehog configuration of the fields in eqs.([S]) also contains isospin degrees of freedom and in principle different 
configurations of solitons in the isospin space can lead to lower energy states in the lattice. Firstly we focus on the 
study of the interaction of the two solitons by changing the relative orientation of the hedgehog quills in the regime 
of long distances, through the use of the product ansatz. This approach permits to describe the soliton-soliton in- 
teraction in the range of long/intermediate distances, namely as long as the two B = 1 solitons do not strongly overlap. 

Assume we have two solitons, each a B = 1 hedgehog solution, whose centers are at ri and r2. For simplicity we 
place the centers of the two hedgehogs symmetrically along the z axis at a distance d, hence the explicit expressions 
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for the centers becomes 

T-n = m n - 

2 



n = (o,o,-^) (21) 



2 

In the product ansatz scheme the new meson fields configurations read: 

u 



r2 = (0,0,^). (22) 



ai{r - n) + ir ■ TVi{r - ri)\ / (J2{r - t-2) + ir ■ TV2{r - | ^^3) 



where /^r is introduced in order to have the proper dimensions and to recover the correct asymptotic behaviour. The 
operator A is a SU (2) matrix inserted to take into account the relative orientation in the isospin space of one of the 
solitons with respect to the other. More explicitly the meson fields for the B = 2 system are given by: 

<JB=2{ri,r2) [(J{ri)a{r2) + 7r(ri) • n{r2)] 

•n-B=2(r"i, ^2) [(j{ri)TT{r2) + a{r2)n{ri) + 7r(ri) x 7r(r2)] (24) 

In the Skyrme model, taking the product of two B = 1 soliton solutions is one of the most convenient ways to obtain 
the B = 2 intersoliton dynamics ^] . The relevant difference between the CDM and Skyrme model is the presence of 
quark degrees of freedom, that also needs to be taken into account in the new fields configuration. At zeroth order, 
since the meson background configuration shows a reflection symmetry except for an isospin rotation, the quark flelds 
can be expressed as a linear combination given by: 

V'±(r) = -i=(Vi(r)±^fl(r)) , (25) 



where the left and right spinors read: 



tpL{r) = ijji{r - ri) , 

^R{r) = AMr-r2). (26) 



In order to study the soliton-soliton interaction as a function of the isospin, we decide to keep one soliton fixed and 
rotate the other, namely the one centred in r2. We will consider three different configurations, shown in Fig. [TJ 

1. Configuration A: A ~1, i.e. the two solitons are unrotated (panel A in Fig. [ij. 

2. Configuration B: A = e*^'^ = ztz, which corresponds to rotating the second soliton by an angle tt about the 
axis parallel to the line joining the two centers (panel B in Fig. [IJ. 

3. Configuration C:^ = e*"^'^ — ir^, leading to a rotation of 180 degrees around the axis perpendicular to the line 
joining the two solitons, in our case the x axis (panel C in Fig. [J). 

In the three configurations A, B and C just one pion is rotated, so we can interpret this transformation as a changing 
in the isospin orientation of the pion. The quantum numbers of the pion are S = 0, I = I and intrinsic parity P = —1, 
leading to a grand-spin = 1~ which admits three different projections Gs ~ —1,0, 1. Hence by rotating the pion 
field we obtain a linear combination of G3 projections. 

The total grand-spin for the B = 2 system is given by the sum of the grand-spin Gi for the soliton placed in ri and 
of G2 for the one centred in r2. Moreover we have to take into account the total parity Ptot of the system. For the 
quark state, since all quarks lie in the 0+ state, they possess positive parity. The pion couples to the spinor through 
the 75 matrix, which under parity changes sign, thus the globally parity carried by the pion coupling will again be 
positive. This leads to a total parity Ptot equal to -1-1. The total grand-spin will follow from: 



"-^2 ~ '^2,quarks "t" '^2,pion ~ ^ 



G^r=0+,l+,2+ (27) 
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A 
B 




Figure 1. The three configurations A, B, and C corresponding to the different orientation in the isospin space. 



In particular, from eqs. (|24|) . it can be seen that the product ansatz represents a possible way to write the pion 
component of the B = 2 solution as a vector which is associated to the quantum number G = 1. The difference 
between the three configurations lies in the different spatial-isospin arrangement of the fields. For instance the only 
configuration that keeps unaltered the hedgehog structure when the two solitons are on top of each other is the A 
configuration (for which Gs^tot — 0). The remaining configurations B and C, due to the rotations, will generate 
a disposition of the chiral fields which will give a different correnspondance between spatial and isospin degrees of 
freedom. 

For each isospin configuration (A,B,C) we calculate the energy of the B = 2 soliton using the approximated field 
expressions given in (j23p and evaluating the expectation value of the Hamiltonian on the new B — 2 state where the 
valence quarks fill the levels provided by eq. ([SS]) . The energy hence reads: 

(28) 



Eb=2 — Ea,B=2 + £^7r,B=2 + Epot,B=2 + Eq^b=2 

where the meson and potential energy are given by: 

Ea,B=2 = ^ I d^r'VaB=2 ■ V(TB=2 , 



Ett 



,B=2 



E.^ 



pot,B=2 



(rrVTTB=2 ■ ^TTB= 



J d^rV{aB=2,TTB=2) , 



B = 2\ d^rilj[ia ■ V + grr{(TB=2 + ir ■ tvb=2)]->P\B ^ 2) 



while the quark energy reads: 

Eq^b=2 

The soliton-soliton interaction energy as a function of the intersoliton separation d is defined as: 

Vint{d) — Eb=2 — 2Eb=1 

where the energy for the B = 1 system is given by eq. ( [T5|) . 



(29) 
(30) 
(31) 

(32) 
(33) 



In Fig. [2] we present the total interaction energy for the three configurations. There is no significant difference in 
shape between the three cases but the three curves show very different behaviours as d decreases. First of all we 
notice that at large and intermediate distances the B — 1 case is recovered, as expected since the two solitons are 
well localized and far apart. As the separation reduces the quark Pauli repulsion starts acting and at d < 1 fm all 
the configurations show a strong repulsive trend. The A configuration is basically repulsive, showing only a small 
attraction of a few MeV around d = 1.5 fm. The B configuration is also fundamentally repulsive except for a small 
attraction at smaller distances, d = 1.2fm. Finally the C configuration is strongly attractive above d — 1.4 fm allowing 
for a strongly deformed bound state. 
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d (fm) 

Figure 2. Soliton-soliton interaction energy as a function of the distance d for the three configurations A, B and C. 



The contribution of the quark fields plays a crucial role at short separation. The large repulsion is due to the fact 
that we are forcing the six quarks to stay in a Grand-Spin = O"*" and does not allow to get a clear insight in 
the short range regime. This area should be explored by using a different approach which allows to deal with the 
fcrmionic degrees of freedom, by including also the excited states for the quarks. 

B. Short distances: the six-quarks bag approach 

In order to take into account the effect of Fermi statistics, we need to include the possibility that quarks fill excited 
states. As the distance between the two solitons reduces, we can imagine that the system would preferably fall into a 
lowest energy configuration in which the six quarks fill the appropriate excited states. At zero separation the B — 2 
system is given by the self-consistent solution of the Lagrangian in eq. (jSj with six quarks. The exact solution of this 
six quarks bag provides the energy of the B = 2 system when d = 0. 

Following [3JJ, we solved self-consistently the field equations for the model with six quarks by assuming the 
hedgehog ansatz for the fields. The hedgehog structure requires a symmetric occupation of the third components 
of the Grand-Spin G; from this constraint the source terms for the a and the pion field are invariant and the field 
equations will differ from the three quarks case just for some multiplying constants. From these considerations, we 
can anticipate that, in order to preserve the hedgehog configuration, the matching between the short and the long 
range interactions will only involve the A configuration of the product ansatz. 

The filling of the quark states requires a brief explanation. The ground state of the hedgehog solution is given by: 

|G^-0+>-|G = G3 = 0,j-i,? = 0,fc = -l) 

where j,l,k are respectively the total angular momentum, the orbital angular momentum and the Dirac quantum 
number. Here we can put three quarks since we include the colour quantum number. 

The next level is the state \G^ = 0^) ^ \G = G3 — 0,j = ^,1 = l,k = 1) and again here we can put the other 
three quarks. The next state in energy available will be the \G^ — 1+) which can accept up to nine quarks since 
G3 = —1, 0, 1. So we will have two possibilities to place six quarks: 

1. First configuration: 0+,0^ with total grand-spin Gf*f = 0~; 

2. Second configuration: 0+, 1+ with total grand-spin G^'j"' = 1+. 



8 



As already shown in ^ijthe set that provides the lowest energy state is the second one, where we place three quarks 
in 1+ ^3] and the others three quarks in 0+. 
The total energy for the six quarks bag reads: 



pot 



(34) 



E = 7VG=oeo + iVG=iei + E„ + E^ ■ 
= 3(eo + ei) + Ed + Et: + Epot 

where eo, ei are respectively the quark eigenvalues for the 0"^ and l"*" states. In Table HIl we show the results for the 
values — 550 MeV and g = 5. we also present the interaction energy Vint — E^q^arks ~ '2Eb=i 



Table II. Contributions to the total energy. All quantities are in MeV. 



Quantity 


Our Model 


eo 


-30.6 


£1 


344.8 


Ea (mass+kin.) 


470.6 


Ett (mass+kin.) 


643.8 


Potential energy a — n 


251 


Total energy 


2308.2 


Total energy B = 1 


1157.8 


Vint 


-7.4 



The next step is to connect the result at d = with the results coming from the product ansatz at large distances. 
We showed that the B — 2 system, written through the product ansatz, admits G^*t°* — 1"^; in addition the solution 
given by the A configuration is the only one that keeps a hedgehog structure. Hence we can expect that as d decreases 
the system will move from the configuration A in the product ansatz to the state provided by the six quarks bag. 

In order to obtain the energy of the six quarks bag as a function of the intersoliton separation, we adopt a variational 
approach, namely we use the exact solutions for the six quarks problem as initial ansatze and evaluate the total energy 
at each value of d. The trial wave functions (following [38]) for the meson fields read: 



cr(r, Rc) = -U cos 



TT tanh 



log3^ 
2i?c 



7r(r,i?c) = -f TV sin 



TT tanh I — ^—r 
2Rc 



(35) 



where the variational parameter Rc is related to the range of action of the attractive potential of the chiral fields. 
The appropriately normalized Dirac components are given by the following expressions: 



u{r, Rd) 



1 



uq exp 



Rjj 



u(r, Rd) 



1 



zVqt exp 



(36) 



VN{Rd) 

Here the parameter Rd is strictly connected to the size of the bag, and hence to the distance d, through the relation: 

d={<rlg >Y'^Rd) = (^j r\u\r,RD)+v\r,RD))dr^ ' 



(37) 



The total energy of the bag depends on the variational parameter R^ Rd and on the quark-chiral fields coupling 
g and is given by: 

E(Rc, Rd, g) — Eint{Rc, Rd, g) + Ekm. quarks {Rd) + Ea{Rc) + E-^{Rc) + Epot{Rc) 

where Eint [Ro, Rd , g) represents the quark-meson interaction energy. As already mentioned previously, the interaction 
energy for the six quarks bag as a function of the distance is defined as: 



V^ntAd) = EiRc, d{RD), g) - 2Eb=i. 



(38) 



The variational method adopted to calculate Vintfi consists in fixing g 
functional in eg . ( [55)1 at each value of Rd- 



5 and varying Rc to minimize the energy 
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IV. RESULTS 



In this section we present the resuhs obtained by studying the interaction energy of the B — 2 system at large 
distances by using the product ansatz and at shorter distances by adopting the description of the six quarks bag. In 
Fig.[2]we present the interaction energy as a function of the distance d between the two sohtons. The behavior of the 
B = 2 system moving from large to small separations, is as follows: 

• at distances d > 1.5 fm, the A configuration, as expected, provides a solution close to summing up the two 
solitons B = 1; at large separation the behaviour of the soliton interaction seems to be similar to that obtained 
in the Skyrme model since when two solitons are far apart, the interaction between them is mediated by the 
chiral fields; 

• as d is reduced the quark repulsion of the product ansatz makes the interaction repulsive until it matches that 
of the six quark bag, thereafter the system undergoes a transition from A to the six quarks bag configuration 
at roughly dAe ~ 1 fm and the system becomes bound. 




^8.0 0.5 1.0 1.5 2.0 2.5 

d(fm) 

Figure 3. Interaction energies as function of the intersoliton distance; here we show the A configuration obtained with the 
product ansatz approach (dashed lines) and the results from the six quarks bag for Rs — Rp — 1.04 (black solid line) 



In the intermediate distances range, around d « 1.5 fm, the A is lightly attractive by a few MeV, but as d gets smaller 
the quark repulsion increases and the interaction energy starts to rise until a crossing between the A configuration 
and the solution of the six quarks bag occurs. This takes place for d ^ 1 fm. For lower values of d the systems goes 
into a stable configuration of six quarks with i?^ ~ 1 fm. From then on the B — 1 hedgehogs are captured to form a 
quasi bound state leading at zero separation to the local minimum given by the self-consistent solution in Sec. IIIIBI 
Thus as the inter-soliton separation shrinks the system falls into a six quarks bag bound state where the inclusion 
of excited states for the quarks has produced an attractive configuration. It should be noted that the optimal value 
of the variational parameter Rc, around 1 fm, is compatible with the usual size of a soliton; moreover, since the 
transition to the six quarks bag occurs at a value of d compatible with Rc, we notice that the overlapping of chiral 
fields does not affect strongly the short range behaviour, which is completely determined by the quarks contribution. 



V. CONCLUSIONS 



We have used a Lagrangian with quark degrees of freedom based on chiral and scale invariance to study the system 
with baryonic number B = 2 and have gone beyond the Wigner-Seitz approximation. We have described how the 
interaction energy behaves as a function of the distance d between the two solitons. When the two solitons are very 
far apart, we have adopted a product ansatz solution. This method allows to take into account the isospin degrees of 
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freedom of the fields and allows an evaluation of the energy of the B = 2 system for different isospin configurations. 
In the long distances regime, the system admits a strongly attractive configuration (C) when one soliton is rotated of 
TT along an axis orthogonal to the line joining the two centers of the B — 1 hedgehogs. The other two configurations 
are very slightly attractive to allow for bound states at large hedgehog separations. In the product ansatz calculation, 
by reducing the inter-soliton distance we are actually forcing six quarks to lie in the ground state with Grand-Spin 
GP = 0+. Since this state allows to accommodate only three quarks as the two solitons approach each other all the 
three configurations develop a repulsive core associated to the overlap of the B = 1 solutions and due to the action 
of the Pauli principle on quarks. To overcome this problem, we have studied the formation of a stable six quarks bag 
at shorter distances. More specifically, we have found the exact solution for the six quarks bag at zero separation 
by solving self-consistently the field equations and by assuming a hedgehog ansatz for the fields. We have found a 
stable six quarks bag bound state at d = with a binding energy of ~ — 8 MeV. Starting from these exact field 
configuration, we have built an energy functional Vmtfiid) and extended the results to finite values of d in order to 
reach the intermediate interaction range { d ~ 1 fm) where the product ansatz description still holds. In order to 
maintain the hedgehog structure of the fields, used in the six quarks bag approach, the only configuration at large 
distances which admits hedgehog structure is the A state, with no relative isospin orientation between the two solitons. 

We have presented the full description of the interaction energy as a function of the inter-soliton separation. The 
B = 2 system, starting at large values of d in configuration A, increases its energy as the distance shortens and 
the repulsion of the quarks becomes dominant. Around d '-^ 1 fm it becomes favorable for the quarks to jump into 
excited states leading to a six quarks bag which provides the most stable configuration and finally the two solitons 
are captured forming a six quark hedgehog state. 

The work presented here can be improved in several ways. The analysis of the B — 2 system, still in this simplified 
version, can be extended to include configurations beyond the hedgehog ansatz, which could lead to lower energy state 
configurations, as already pointed out in previous works [33]. The study of the B = 2 system and more in general of 
a multi-soliton can be performed by building a three dimensional lattice with specific symmetries that leads to lower 
energy configurations [28|. Moreover, since at the moment we are performing calculation at zero temperature, the 
next step would be to include the dynamics of the dilaton in the model and perhaps also include the vector mesons 
in the calculation in order to perform the study at finite density and temperature, as in Ref . [1^ [13] ■ 
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